Multiscaling, Ergodicity and Localization in Quasiperiodic Chains 
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We report results of numerical simulations of wave-packet dynamics in a class of chains consisting 
of two types of weakly coupled clusters arranged in a quasiperiodic sequence. Properties of eigen- 
states are investigated using perturbation theory of degenerate levels in the coupling strength v and 
by numerical diagonalization. Results show that wave packets anomalously diffuse via a two-step 
process of rapid and slow expansions, which persist for any t; > 0. An elementary analysis of the 
degenerate perturbation expansion reveals that non-localized states may appear only in a sufficiently 
high order of perturbation theory, which is simply related to the combinatorial properties of the 
sequences. Numerical diagonalization furthermore shows that eigenstates ergodically spread across 
the entire chain for w > 0, while in the limit as u — » ergodicity is broken and eigenstates spread 
only across clusters of the same type, in contradistinction with trivial localization &t v = Q. An 
investigation of the effects of a single-site perturbation on wave-packet dynamics shows that, by 
changing the position or strength of such an impurity, it is possible to control the long-time wave- 
packet dynamics. By adding a single impurity it is possible to induce wave-packet localization on 
individual subchains as well as on the whole chain. 
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Understanding the relation between spectral proper- 
ties of a given Hamiltonian and the dynamics of wave- 
packets governed by it remains one of the elementary 
questions of quantum mechanics that still poses signifi- 
cant challenges, further emphasized by the discovery of 
quasicrystals P, 01 . While spectra of many Hamiltonians 
decompose into a point part and an absolute continuous 
part with, respectively, bound (localized) and unbound 
(delocalized or extended) eigenstates, there is a large va- 
riety of Hamiltonians whose spectra, for certain values 
of the parameters, are neither pure point nor absolute 
continuous nor a combination of both. In this case the 
spectrum contains a singular continuous part, with multi- 
fractal eigenstates and/or Cantor-set spectra. Examples 
are Harper's model of an electron in magnetic field 0, 
the kicked rotator |^, as well as the Anderson model 
of an electron in disordered medium 0]. In the case of 
an electron in a quasiperiodic system, as studied in this 
Letter, many examples lead to spectra which are purely 
singular continuous [l(Tj |. 

It is known that a particle's return probability decays 
with the same powerlaw as the scaling of the local density 
of states 10, and that the spreading of the wave-packet 
width d{t) (defined below) exhibits anomalous diffusion, 
d{t) ^ t^, with 0</3<l(/3 = corresponds to the 
absence of diffusion, /3 = 1/2 to classical diffusion and 
/3 = 1 to ballistic spreading) In addition, the wave- 
packet dynamics exhibits multiscaling, where different 
moments of the wave- packet scale with different, non- 
trivially related exponents /3 H Q E HI 113 ■ While 
wave-packet localization implies a pure-point spectrum, 
the converse is not true, and the more refined notion of 
semi- uniform localization is necessary |2nj . The exact 
relations between particle dynamics and singular or ab- 
solute continuous spectra, on the other hand, is less well 
understood. As a rule, systems with singular continuous 



spectra exhibit anomalous diffusion, while absolute con- 
tinuous spectra may lead to either anomalously diffusive 
or ballistic dynamics pj] |. 

In this Letter, we study eigenstate and wave-packet 
properties of an electron moving in a one-dimensional 
quasiperiodic system described by the tight-binding 
Hamiltonian H = + H.c). Here Cj(c]) 

is the annihilation (creation) operator of an electron at 
site j, and the nearest- neighbors-only hopping integrals 
tjji take values 1 or < v < 1 for, respectively, let- 
ters / or s ( "large" and "small" hoppings) of a quasiperi- 
odic sequence of letters obtained by the inflation rules 
{s |—^ 1,1 i—> lsl"~^} iterated m times, starting with the 
letter s. Denoting the length of the m-th iterant by Nm, 
these sequences have the property that Nm/Nm+i 5n 
as m — > oo, where <5„ is an irrational number with con- 
tinuous fraction representation [n, n, n, . . .]. The quan- 
tities of interest are the solutions of the eigenproblem 
HijJk{j) = Ekipkii), k = l,...Nra + l, as wefl as the 
width d{t) of a wave-packet '^{j,t), 
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which is initially localized at a site jo, so ^(j, t = 0) = 

Figure 1 shows the change of the wave-packet width 
with time for the silver-mean (n = 2) model. While 
the asymptotic behavior can be characterized as anoma- 
lous diffusion d{t) , f3 ~ 0.2, there are also "flat" 
parts in the regime of strong quasiperiodic modulation 
(small v). As demonstrated by the insets in Fig. 1, 
these can be characterized by the existence of time inter- 
vals, exponentially growing at a constant rate for a given 
V, during which d(t) strongly oscillates in a self-similar 
manner reflecting the hierarchical structure of system. 
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FIG. 1; Evolution of the width d{t) of a wave-packet initially 
localized in the middle of the sliver-mean chain for m — 9, 
for several small values of v. The insets show two magnified 
"steps" for 11 = 0.025. 

but nevertheless remains bounded from above by a con- 
stant. Further inspection of the wave-packet behavior 
reveals that breathing modes are responsible for the oscil- 
lations, while the wave-packet spreading itself is limited 
to low-amplitude "leaking" out of the region in which it 
is confined. Eventually, the wave-packet expands fast to 
reach the next level (this spreading can be described as 
d{t) (X with P' « 0.79), and the whole process repeats. 

Such a behavior is in an agreement with results of 
Ref. 16], where, based on a qualitative model of the 
wave-packet spreading in the seniiclassical approxima- 
tion and numerical simulations, it was argued that an 
hierarchical splitting of the spectrum into constant-width 
bands leads to a step-like behavior with /3' = 1, which 
then is smoothed out due to the (broad) distribution of 
band widths. The here observed value /?' < 1 for small v 
suggests that there is a dispersion of band widths even in 
the hmit as ?; ^ (when band widths go to 0|), and 
that the self-similar spreading of the wave-packet is only 
an approximate description of a more general multiscale 
dynamics. This self-similarity can be used, e.g., for cal- 
culating P hy (3 ~ Slnd/Slnt, where d\nd and Slut are 
the horizontal and vertical displacements between each 
of the two steps in Fig. 1, respectively, giving in the case 
studied /3 « 0.20 for v = 0.025. 

The self-similarity of quasiperiodic sequences was pre- 
viously used in a renormalization-group perturbative ex- 
pansion that provided a great deal of insight into the 
eigenstate properties 17|, and showed multiscaling of 
wave-packet dynamics :21;1. Here we focus on ergodic 
rather than hierarchical 18] properties. By an elemen- 
tary analysis of the perturbation theory of degenerate 
levels at u = for small v we show that (a) for any w > 
eigenstates delocalize, in contradistinction to (trivial) lo- 



calization at V = 0; and (b) in the limit as u ^ eigen- 
states delocalize across only one set of clusters containing 
the same number of atoms (subcluster localization due to 
ergodicity breaking). 

Raleigh-Schrodinger theory allows the recursive con- 
struction of matrices in subspaces of a given degenerate 
eigenenergy to a given order p, whose diagonalizations 
(the "secular problem") yield corrections to the unper- 
turbed eigenenergies, hopefully accurate up to 0{vP~^^). 
In the case of an unperturbed chain {tjji = t), for in- 
stance, the solution of the secular problem in first order 
in t is equal to the exact solution of the problem. 

For the quasiperiodic sequences considered here, how- 
ever, such perturbation expansions yield two qualita- 
tively different types of solutions, depending on the val- 
ues of p and n. To that end we first notice that ap- 
proximants for a given n consist of two types of words, 
l"^^ and Z", separated by the letter s, corresponding to 
clusters with n and n + 1 atoms, respectively, coupled 
via hopping of strength v which is treated as a pertur- 
bation. The unperturbed system then has 2n -|- 3 degen- 
erate levels with all eigenstates localized on individual 
clusters. In higher orders of perturbation theory, these 
localized states spread as the coupling among the clus- 
ters of the same type is taken into consideration, and, 
for a sufficiently high order, delocalize across the whole 
chain. Since the maximal number of letters s between 
two consecutive clusters of length k is also fc, the eigen- 
states of the clusters and delocalize only in the n 
and n + 1 order of the perturbation theory, respectively. 
More precisely, the dimension of the secular problem for 
each cluster of length q = n,n + \ changes from at most 
O(n^) for p < qto 0{N,,n) for p > q. Only the latter case 
allows for multifractal and/or extended states. 

As an example we analyze the case n = 2. There 
are 7 unperturbed levels, with E^^ = ±(V5 ± l)/2 and 
Eii^ = ±\/2, 0. In first order, the only correction are 
6 bands linear in v, splitting off the three E^"^ bands, 
because sllsUs is the only possibility where two clusters 
of the same type are connected by a single s bond. For 
p ~ 2, all U clusters become connected while III states 
are still not extended due to the existence of lllsllsllslll 
sequences in the chain. For p = 3, states belonging to III 
bands delocalize as well. 

To investigate further the issue of convergence, we no- 
tice that, even when calculated to all orders of v, the 
secular problems for the two types of clusters give solu- 
tions that are inevitably localized on the clusters of the 
given type and with zero component on the cluster of the 
other type. For various values of we check by numerical 
diagonalization whether this is confirmed. Fig. 2 shows 
Pui{Ek) = the total probability that the 

particle in an eigenstate with an energy Ek will be 
on the cluster III. The result implies that pm strongly 
depends on the energy band of the given eigenstate for 
small V, being large for the states belonging to the III 
bands, and vanishing for the states from II bands in the 
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FIG. 2: Probability that a particle in an eigenstate with the given energy is on the III sublattice, for v = 0.1, 0.2, . . . 0.9 (going 
from top to bottom in the central part of the plot). The line connecting points for v = 0.5 is for guidance only. Only states 
E < are shown due to the symmetry about the band center. 



limit V ^ 0. While Fig. 2 shows results for m = 9, we 
have also checked pm for the iterants m = 10, 11, 12 to 
see whether there is any systematic deviation from the 
values shown in Fig. 2. We found that, while there are 
many more states for larger m's, the values of pm do not 
shift systematically for almost all of the states . Instead, 
the additional points cluster in the same way as in Fig. 2. 

This supports the possibility that eigenstates remain 
ergodic (i.e., spread over both types of clusters) even for 
very small values of v, and that, in turn, the Schrodinger- 
Raleigh perturbation expansion of degenerate levels has 
zero radius of convergence but, nevertheless, might still 
be accurate in the limit as u — > since pm (u — > 0) ^ or 
1 (which is, however, only a necessary but not a sufficient 
condition). On the other hand, the nearly non-ergodic 
behavior of wave-packets for small v has interesting con- 
sequences when a single impurity is added, as we show 
next. 

An impurity H' ~ Uic\ci added at a site i will, for 
large u {e.g., it ^ 1), act as a barrier, effectively cutting 
the chain into two halves. The consequence of ergodicity 
in this limit is that the wave-packet reflects off the impu- 
rity independently of its initial site, even if v is small. For 
M — > 0, on the other hand, the unperturbed wave-packet 
propagation of the case u = is restored. Understanding 
the wave-packet propagation in the regime of interme- 
diate u values, however, still poses significant challenges 
and surprising results [20| ■ We performed several numer- 
ical experiments for various initial positions of impurities 
and wave-packets for intermediate values of u, and one 
common situation is shown in Fig. 3. The evolution of an 
initially localized wave-packet exhibits high sensitivity on 
the position of the impurity, approaching two quite dif- 
ferent stationary states. In particular, while in Fig. 3(a) 
the final state is just slightly perturbed from a final state 
for u = 0, in Fig. 3(b) certain parts of the wave packet 
get reflected, with a small-amplitude "leak" reminiscent 
of the similar process mentioned above. This kind of 



wave-packet dynamics is a consequence of the nearly non- 
ergodic properties of eigenstates discussed above (cp. Fig- 
ure 2). 

Finally, we address the influence of an impurity at the 
initial site of the wave-packet on its dynamics, by study- 
ing the dependence of the flnal wave-packet width on the 
value of u. Fig. 4 shows the maximum value of d{t) at- 
tained in the course of the evolution of a wave-packet, 
which was initially localized at the site = Nm/'^., in 
the presence of an impurity located at the same site. The 
m — 7 iterant of the silver-mean model has unperturbed 
wave-packets exhibiting one full step from Fig. 1, reach- 
ing the final stationary state after beginning of the second 
step. 

Fig. 4 shows that, apart from the strongly localized 
final wave-packets for large u and unperturbed final sta- 
tionary wave-packets for small u, there is a wide range of 
values of u for which the final width of the wave-packet is 
significantly reduced even for u <ti v, signaling (dynami- 
cal) localization. There are nevertheless several well de- 
fined peaks about some values of u where the maximum 
wave-packet width is significant enhanced compared with 
the case for slightly smaller and slightly larger values of u. 
We additionally checked whether such peaks persist when 
the system size is increased, by performing the same cal- 
culation for the TO = 9 iterant. Fig. 4 shows that there is 
again a peak structure, with both small and large peaks. 
The relation of the values of u where these peaks ap- 
pear and the spectrum of H remains unclear. Repeat- 
ing similar numerical experiments for various values of u, 
we find that the peak structure persists as long as the 
band widths are smaller than the band gaps (v ^ 0.4 
in the silver- mean case). These results, however, show 
that in quasiperiodically modulated quantum wires one 
can strongly influence the electronic transport by induc- 
ing local perturbations that act as sort of control gates 
for long-range quantum interference effects in electronic 
transport. 
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FIG. 3: Snapshots of the evolution of two wave-packets, both starting as (5-functions in a local environment . . . sl{x)lls ... at 
the same initial site (a;), in the presence of an impurity. In the two panels (a) and (b), the impurity is located in a local 
environment . . . sll{b)sll{a)slls ... at the position marked by (a) and (b), respectively, as indicated in each panel by a vertical 
line. For easier comparison, the vertical dashed line in each panel marks the position of the impurity in the other panel. The 

long-time wave-packet dynamics exhibit high sensitivity on whether the impurity is located or not on the cluster of the same 
type as the one where the wave-packet was localized initially. The remaining parameters are n = 2, m = 8, v = u = 0.1. 
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FIG. 4: Maximal width of a wave-packet attained during its 
evolution in the presence of a single impurity of strength u 
located at the initial site of the wave-packet. 
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